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We propose a simple variational form of the wave function to describe the ground state and 
vortex states of a system of weakly interacting Bose gas in an anisotropic trap. The proposed wave 
function is valid for a wide range of the particle numbers in the trap. It also works well in the 
case of attractive interaction between the atoms. Further, it provides an easy and fast method to 
calculate the physical quantities of interest. The results compare very well with those obtained by 
purely numerical techniques. Using our wave function we have been able to verify, for the first time, 
the predicted behaviour of aspect ratio. 

PACS numbers: 03.75.Fi, 03.65.Db, 5.30.Jp 



I. INTRODUCTION 

Observation of Bose- Einstein condensation in cooled and trapped dilute gases of alkali atoms jl|-f| and spin polarized 
atomic hydrogen Q has generated a renewed theoretical interest in understanding such systems. In a meanficld 
approach, which is valid in the limit pa 3 << 1, where p is density of atoms and a is the s-wave scattering length, 
ground state and vortex states of these systems can be described by Gross-Pitaevskii (GP) equation [^J. Various 
numerical procedures |6|-[l0[ and approximate analytical method s |pd| - [l5[| have been used to solve the GP equation. 
Among these variational scheme proposed by Baym and Pethick Jll| to explain the experimental observations of Ref. 
jl| is particularly appealing. In this approach the trial wave function was taken to be of the form of ground state of 
the trap potential (modeled by an anisotropic harmonic oscillator potential). Thus the wave function is represented 
by a three dimensional Gaussian with axial and transverse frequencies as variational parameters. This form of wave 
function, however, is valid only when the number of atoms in the trap is very small. As the number increases, the 
repulsive interaction between the atoms tends to expand the condensate and flatten the density profile in the central 
region of the trap where the density is maximum. Of these two effects, only expansion of the condensate can be 
described adequately by the Gaussian trial wave function. On the other hand, the Thomas-Fermi approximation 
IPIptI provides a wave function which is valid when the number of atoms is very large which is the case with the recent 
experiments ]l6| . It is of interest, however, to have a form of wave function which apart from showing these limiting 
behaviours, is also valid in the intermediate region. In a variational approach, Fetter |l4| has proposed such a wave 
function for the bosons in an isotropic trap. 

In this paper we propose a simple form of wave function for the ground state of bosons confined in an anisotropic 
trap. The wave function is valid for a very wide range of particle numbers. When the number is small it tends to 
mimic a Gaussian, and in the opposite limit it resembles the Thomas-Fermi wave function. However for large number 
of atoms the wave function differs from the Thomas-Fermi wave function in the surface region, a desirable feature as 
noted in Ref. |9|,[l7j . The trial wave function has got an additional parameter comapared to the ones used in Ref. [jll| . 
This parameter takes care of the flattening of the density in the central region of the condensate. Thus, it provides a 
better lower bound for the ground state energy than the Gaussian trial wave function. We also compare the results 
obtained by our trial wave function with those obtained by other numerical procedures such as the minimization of 
energy functional by steepest descent method |J and the integration of the nonlinear Schrodinger equation HpXfl . 
These comparisons show the form chosen by us to be highly accurate for obtaining a host of physical qauantitics 
of interest. In addition to providing acurate results for a very wide range of the particle numbers, the method is 
also very fast from the computational point of view. Further, using this wave function, the physical observables 
can be expressed analytically in terms of three variational parameters which are obtained by minimizing the GP 
energy functional. Since we have a simple analytical form for the energy in terms of three variational parameters the 
procedure of minimization is very simple. The novel achievement of this method lies in verfication of the predicted 
behaviour of the aspect ratio which is very important quantity from the experimental point of view. It could not be 
ascertained before because convergence of the aspect ratio to the highly repulsive limit is very slow Since our 
method can handle even very large number of atoms in the trap without any difficulty we could verify the behaviour 
of the aspect ratio. 
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Based on physical reasoning we generalise the variational form to descibe the vortex states also. As is the case with 
the ground state, we find good agreement with the existing results, with considerably less computational effort. 

Condensation has also been observed in 7 Li [Q where interatomic interaction is attractive which is characterized 
by the negative s-wave scattering length. As the number of atoms in the trap inceases the condensate shrinks and 
nonuniformity in the central region increases. After a critical number of atoms in the trap, the condensate collapses. 
This situation is also very well descibed by our trial wave function. The parameter which accounts for flattening of 
the density profile in the case of repulsive interactions also takes care of the increase in density gradient in an effective 
way . The critical number of atoms for the case of 7 Li is in close agreement with that reported in Ref . |fio|| . 

The paper is organized as follows. Section II contains the description of the variational scheme employed in the 
paper. It also contains the analytic expressions for the observables of interest. Results obtained from the variational 
procedure and their comparison with the existing ones are presented in Section III. Section IV contains summary and 
concluding remarks. 



II. VARIATIONAL METHOD 



A. Ground state 



Bose Einstein condensation in experiments with cooled and trapped atoms can be described within the framework 
of the GP theory. Validity of such a description has been analysed by Stenholm Jl8| . In a situation where the trap can 
be modeled by an anisotropic harmonic oscillator potential characterized by the two angular frequencies to^_ and 
and the interatomic ineractions can be replaced by an effective pseudo-potential involving s-wave scattering length a, 
the ground state energy for condensed bosons of mass m is given by the Gross-Pitaevskii functional || 
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|V^i(n)| 2 + {x\ + y\ + Xlzl) Wn)! 2 + f |V>i(ri)| 4 j . (1) 



Here we have used the length scale a±_ = ^h/muj^ and the energy scale Kuj^ provided by the trap potential to 
express Eq.(Q) in terms of the dimensionless variables. V'i( r i) is the condensate wave function which satisfies the 
normalization condition 

'drihMrOl^l. (2) 

Ao = u>^/u>j_ is the anisotropy parameter of the trap, N is the total number of atoms in the condensate and u± = 
8iraN/a±. The exact form of the wave function can be determined by minimizing the energy functional in Eq.(Q) 
with the normalization constraint of Eq.(Q). Such a minimization results in the nonlinear Schrodinger equation 

[-V 2 + (xl +v\ + Xlzl) + Ul |Vi(rO| 2 ] Vi(ri) = 2 Ml Vi(r 1 ), (3) 

where, [i\ is the chemical potential. It is not possible to find exact analytic solution to Eq.(||). Consequently various 
numerical techniques have been developed to study the ground state property of such systems within the framework 
of the GP theory. These techniques involve either the direct numerical minimization of Eq.([l]) with the constraint 
of Eq.(|J) U or numerical integration of Eq.(Q) or its time dependent version ]q-PJTo|. Another approach is to use 
the variational method which has been extensively used in different branches of physics. The main advantage of this 
method is that with a suitable guess for the form of the wave function it is possible to save a lot of computational 
effort and time. In addition, it may also provide physical insights which generally get obscured in the complicated 
computational procedures. The first study of this kind was done by Baym and Pethick in light of the experimental 
observations in 87 Rb jij. They took the trial wave function for the ground state as 

V,(r) = ATV2 1/3^1/4 f ™) e -m(u,xrl+«.» a )/2» (4) 

with effective frequencies, u>± and lo z , treated as variational parameters. However, the wave function above brings 
out only the qualitative features of the condensate e.g. expansion of the condensate in different directions, shifts 
in the angular frequencies and the scaling behaviour of energy with the number of atoms in the trap. Further, this 
form of the wave function is valid only for small number of atoms in the trap (see Fig. 2 below). We now propose a 
variational form of the wave function and demonstrate its applicability and utility in providing accurate description 
of the condensate for a wide range of the particle numbers. The form of the trial wave function we choose is 
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where, X,to± and p are the variational parameters which are obtained by minimizing the energy E\ in Eq.(Q) with 
respect to these parameters. It is easily verified that the wave function satisfies the normalization condition of Eq.(^|). 
The expression of the ground state energy E\ in terms of A, u>± and p is 

El - T2^{ 1+ 2)-^ { 1 + 2P)+ 3^{ 1+ 2X 



a / to 



V\p (i 



2p 

For a particular value of N the parameters tu± , A and p are obtained by minimizing the energy above using standard 
numerical routines. We have used Mathematica JliJ for this and it takes a few seconds of the real time on 166 MHz 
Pentium-1 computer to get the answer. Next we discuss how different physical observables can be obtained in terms 
of the parameters of the wave function. Aspect ratio which characterizes the anisotropy of the velocity distribution 
of the condensate is defined as y/ < pi > / < p 2 >. This can be easily shown to be 



<pl> = <j4> 

V < pi > V < z i > 

Width of the condensate in the transervse direction is given by 



v/A. (7) 



and the width of the the momentum distribution in this direction is given by 



2 < r {£) 

< x\ >= y-4- (8) 



NhmuJiT (i) (1 + 2p) 
<Pl>= ■ (9) 
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The peak density of the condensate is Nipf (0) / a 3 ± . Life time of the condensate is related to the density distribution. 
The loss rate due to the two body loss rate M and the three body loss rate |M[ is given by 



(10) 



R(N)=a J dr^O) | 4 +L J dr\4>(r)\ 6 

_ aN 2 V\(LU ± /Lu a ± ) 3 / 2 p LN 3 X{lo ± /lu° ± ) 3 p 2 
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where, a is the two-body dipolar relaxation loss rate coefficient and L is the three-body recombination loss rate 
coefficient. 



B. Vortex States 



We consider here the states having a vortex line along the z axis. Wave function of such a state can be written as 

#(r) = ?A(r)e IK * (11) 

where k is an integer denoting the quantum of circulation. Subsituting the complex wave function ^ in place of ip in 
Eq. (m we get Gross-Pitaevskii functional for the vortex states in terms of the scaled variables 
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NiMn)\ 2 + (^l + r 2 ± + \ 2 z 2 ) |^( ri )| 2 + ^(rOI 4 



(12) 
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The corressponding nonlinear Schrodinger equation is 

[-V? + K 2 r^l + rf ± + \%zl + ui|^i(n)| 2 ] Vi(ri) = 2/xn/>i(n). (13) 

We assume the trial wave function of the form 

^(n) = Ar\ L e-^^<y^ +Xz *Y (14) 

where q is an additional variational parameter. This particular form of the wave function is motivated by the following 
considerations. 

1. Presence of the centrifugal term K 2 /rf ± forces the wave function to vanish on the z axis. 

2. It has been shown that for a weakly interacting Bose gas [ p2| the wave function corresponding to fcth quantum 
circulation behaves as 

V- ~ r* ± (15) 

near the z axis. 

Proportionality factor in Eq.(|l4|) is determined by the normalization condition (Eq.^J) 

A , Mil!) ^a (3W 

7r3/2r(i + g )r(^) \<A) 

For a vortex line descibed by the wave function in Eq. ([l4]) the density peaks at 




r ^ = \—[ 1 > ( 17 ) 

V ^ \P/ 

and the peak density is given by 

It is also sraightforward to get the analytic expression for E\ in terms of the variational parameters w, A, p and q 
which, in turn, are obtained by minimization of E\. The kinetic energy is given by 

+ 2g) [(1 + 2p)(l + A/2) + 9 (2p + 2g + A)] T (l±*a) 
4(3 + 2?Kr(f) 

Energy corresponding to the rotational motion is given by 

^i/iV) r o t = , V / 7 - (20) 

The oscillator energy is given by 
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(Ei/N) HO = — / Qin .N ■ ( 21 ) 

The interaction energy is given by 



3/2 a VAiVp(i + 2 (Z ) 2 r(2 (Z )r 2 (i + g )r(3±i2) 



{E-i/N) vot =21^) V 2P r . (22) 

2(3+49)/2 P7r i/2 a±g (i + 4 g )r2( g )r(i + 2q)r 2 (2±|£j 



It is easy to verify that ground state is obtained by setting k and q equal to zero. Once we have the energy of the 
states with and without vortices we can calculate the critical angular velocity for the formation of the vortex state. 
In the unit of w° it is given by p 

il c = K- 1 [(E 1 /N) K -(E 1 /N) ]. (23) 

To demonstrate the applicability of this method we have performed calculations for Rb and 7 Li. The s wave 
scattering length is positive for 87 Rb. It is negative for 7 Li. Consequently the interatomic interaction is repulsive in 
the former and attractive in the latter. We now present the results and their comparison with the existing calculations. 



III. RESULTS 



A. Positive scattering length: Rb 

In this section we report calculations on 87 Rb. We perform our calculations by employing the experimental numbers 
for the asymmetry parameter of the trap, the axial frequency and the s wave scattering length corresponding to the 
experimental situation of Ref. Q| and the subsequent theoretical calculations 0^,0. Accordingly, A — ll>®/lu^ = \/8; 
ujz/2ir is 220 Hz; a is 100a , where ao is the Bohr radius. The corresponding characteristic length is ai = 1.222 x 1CP 4 
cm and the ratio a/a± is 4.33 x 10~ 4 . 

First, we obtain the energy E\ (Eq.^|) by minimizing it with respect to the variational parameters u>±, A and p for 
various values of the particle number N. We present the results in Fig.l. It is evident that the results obtained by 
us are in close agreement with the results in Ref. ||(see also Table 1 below). As pointed out above, these agree well 
with the results obtained by using the Gaussian trial wave function when N is small and with those obtained by using 
the Thomas-Fermi approximation ]9|,pj| when N is very large. These comparisons clearly establish the validity of our 
wave function for a very wide range of the particle numbers. 

Next, we compare the proposed wave function with the Gaussian trial wave function of Eq.(^) and also the one 
given by the Thomas-Fermi approximation in Fig. 2. It is clear that when the number of atoms in the trap is small the 
proposed wavefunction has more resemblence with the Gaussian wavefunction (Fig. 2a). As N is increased the wave 
function tends to flatten in the central region and the resulting form is a mixture of the two limiting forms i.e. the 
Gaussian and the Thomas-Fermi wave function. In the central region it is close to the latter while it resembles the 
former away from the trap centre (Fig. 2b). When N is very large the resemblence is more with the Thomas-Fermi 
wave function (Fig. 2c). However we note that the wave function vanishes smoothly far away from the centre of 
the trap. As mentioned above this a desirable feature which is crucial for the calculation of some relevent physical 
observables e.g. the aspect ratio j9|. It is clear that our wave function not only provides a better lower bound for the 
energy but also shows the correct and the desirable limiting behaviour. 

Results of calculation of various quantities e.g. chemical potential, total energy , kinetic energy, potential energy, 
interaction energy, average size of the condensate in the tranverse direction and in the axial direction have been 
presented in Table 1. The close agreement with the results of Ref. is evident. We have also looked at the variation 
of the peak density and the total loss rate of the atoms with the particle number N. We find them to be in good 
agreement with the result of Ref. |s| . 

We now present the behaviour of the aspect ratio which is a very important quantity from the experimental point 
of view . As mentioned in Ref. Pilj| it is equal to v^o m the non-interacting limit and tends to Ao in the highly 
re pul sive limit, which is the case when N is very large. However the convergence to the repulsive limit is very slow 
PJil"| and therefore this behaviour has not been seen explicitly so far. On the other hand, with a variational wave 
function, calculations can be performed for any N with equal ease. Consequently we have been able to verify the 
predicted behaviour of the aspect ratio. It is seen from Fig. 3 that it requires calculations up to N ~ 10 6 to see the 
aforesaid behaviour. 

In Fig. 4 we show the wave function of the vortex state corresponding to n — 1 for N — 5000. The atoms are pushed 
away from the z axis. Peak density is 7.155x 10 13 which occurs at ri± = 1.611. The position of the peak moves further 
away from the z axis as N is increased. For N = 10000 the peak occurs at r%± = 1.844 while it is at ri± = 2.161 for 
N = 20000. Value of the peak density remains much the same for N = 10000 and 20000( 9.524xl0 13 and 12.26xl0 13 , 
respectively). We also compare the energy of the vortex state with different values of k with that obtained in the 
Thomas- Fermi approximation p3| . For calculation in the Thomas-Fermi approximation we have neglected the change 
in the chemical potential due to the presence of a vortex line. Results are presented in Table 2. The Thomas- Fermi 
approximation yields lower value of the energy which is expected as it ignores the kinetic energy contribution from 
the core and the surface regions of the vortex state. This approximation is valid only for very large N. We do find the 
results of the two calculations to be in close agreement for large N. For N = 10 6 the two results differ only by ~ 1%. 
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We have also calculated the critical angular velocity for the vortex state with k = 1 for various values of the particle 
number N. We find that it decreases rapidly in the begining (up to TV ~ 2000). Thereafter it varies slowly. The 
result shows good quantitative agreement with those in the Ref. ^,^|. For example for N — 2000 the critical angular 
velocity is 52% of the noninteracting value given by the transverse angular frequency u>^_ of the trap, in comparison 
with 49.33% in Ref. [pi. For N > 5000 it is less than 43% of the noninteracting value, which compares well with the 
figure of 40% in Ref. The critical angular velocity increases with k. For N = 10000 we find Sl c /2it equal to 30.57, 
38.5, and 45.42 Hz for k =1, 2, and 3 respectively. These figures are 26, 35, and 41 Hz respectively in Ref. ||. 

B. Negative scattering length: 7 Li 

In this section we report calculations on Li. These atoms interact via an attractive interaction and consequently 
the scattering length in this case is a negative quantity. Numerical values of the parameters used in the calculations 
correspond to the experimental situation of Ref. M and the subsequent theoretical calculations |]|l(| . Accordingly, 
the asymmetry parameter of the trap is A = w°/a7[ = 0.72. The axial frequency w°/27r is taken to be 117 Hz. The 
s-wave scattering length a is — 27ao . The corresponding characteristic length is a± = 2.972 x 10 -4 cm and the ratio 
a/a ± is -4.33 x 10~ 4 . 

We find the value of the critcal number N c = 1270 beyond which the ground state collapses because of the attractive 
interaction. This is in good agreement with the figure of N c ~ 1300 in Ref. (To) and the experimental observation (see 
the second paper of Ref. ||) that N < 1300. Wave functions for N = 500 and 1270 are shown in Fig. 5. For N = 500 
there is hardly any difference between the proposed wave function and the Gaussian trial wave function. However, the 
difference is rather significant for N — 1270 as is evident from the figure. We plot the aspect ratio for various values 
of N < N c in Fig. 6. At N ~ N c the aspect ratio tends to 1. Since for a wave function of the form given by Eq.d3) the 
aspect ratio also gives ratio of spatial widths in the transverse and the axial directions, the condensate tends to be 
isotropic for N ~ N c . This becomes further evident in Table 3 where we have listed the results for various quantities 
of interest. As reported in Ref. H the variation in the various quantities is smooth from N = 1 to 1000. However, 
we can also note the sharp variation as we reach the critical number. This behaviour is consistent with that reported 
in Ref. We have also calculated the peak density and loss rate for N < N c . Once again we find a sharp increase 
near N ~ N c . These results also match well with those of Ref. |l(]] . 

It is possible to have very large number of particles (N » 1300) in the vortex states even when the interatomic 
interaction is attractive. We have considered the vortex states with k = 1, 2 and 3. The particle number is 3500, 
6000 and 8000 respectively. Peak densities for these states are 1.266 x 10 13 , 2.239 x 10 13 and 2.744 x 10 13 which occur 
at ri± = 0.922, 1.257 and 1.571 respectively. Although the particle number is quite different in the three cases, the 
peak densities are not very different. Also they remain less than the peak density 3.984 x 10 13 which corressponds 
to k — and N — 1270. These observations are consistent with those in Ref. Q. Stability of the vortex states for 
attractive interaction can be physically explained as the interplay between the restoring force and and the centrifugal 
force. The restoring force tries to attract the particles to the centre while the centrifugal force tries to push them out. 
The net effect is that the peak density does not change much even when there is siginificant variation in the particle 
number. Since the interparticle interaction depends on the density, for low dnesities it does not cause the collapse of 
the condensate. 

In case of attractive interaction it takes more energy to create a vortex state than that requiread in the noninteracting 
case. Consequenty the critical angular velocity is greater than unity. For k = 1 and N = 1000 we find £l c = 1.119 
which compares very well with f2 c = 1.2 reported in Ref. 

IV. CONCLUSION 

We have proposed a variational scheme to describe the ground state and vortex states of weakly interacting atomic 
gases confined by harmonic traps within the framework of the meanfield theory of Gross and Pitaevskii . It is based 
on a judicious choice of the form of trial wave function for the ground state which has a simple functional form and 
at the same time is valid for a wide range of the particle numbers. When the number is small it tends towards a 
Gaussian and in the opposite limit it resembles the Thomas-Fermi wave function. However, for large N it provides a 
better description of the surface region than the Thomas- Fermi wave function. In the intermediate regime it combines 
the feature of both in an effective way. In the central region of the trap, where density is high, it matches with the 
Thomas-Fermi wave function. Away from the centre of the trap, where density is low, it matches the Gaussian trial 
wave function. We easily generalize the wave function for the vortex states. We have demonstrated the applicability 
of our method by performing calculations of various physical quantities for the experimental situations of Ref. fflj and 
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Ref. g. We find our results to be in good agreement with the existing results. The method is semi-analytic and 
consequently computationally easy to implement. As our method poses no additional problems even for very large 
particle numbers we have been able to verify, for the first time, the predicted behaviour of the aspect ratio. The 
formalism is quite general involving only the scaled s-wave scattering length and the asymmetry parameter of the 
trap. We therefore believe that it will be useful in analysing a variety of experiments. In addition it may serve as a 
very good starting point for the theories [p4| where quantum fluctuations play an important role. Generalization of 
our method to the time dependent case ]12| is straightforward and it will be reported in a future publication. 
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TABLE I. Results for the ground state of 87 Rb atoms confined in an anisotropic harmonic trap with Ao = %/8 and 
Wx/27r = 220/Ao Hz. Chemical potential and energy are in units of %uj]_ and length is in units a±. Numbers in the brackets 
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0.755 


2.16 


0.936 


1.078 


0.499 




(4.77) 


(3.84) 


(0.76) 


(2.15) 


(0.93) 


(1.08) 


(0.5) 


2000 


5.951 


4.628 


0.66 


2.645 


1.323 


1.227 


0.534 




(5.93) 


(4.61) 


(0.66) 


(2.64) 


(1.32) 


(1.23) 


(0.53) 


5000 


8.164 


6.142 


0.543 


3.577 


2.022 


1.469 


0.596 




(8.14) 


(6.12) 


(0.54) 


(3.57) 


(2.02) 


(1.47) 


(0.59) 


10000 


10.527 


7.783 


0.461 


4.577 


2.744 


1.689 


0.657 




(10.5) 


(7.76) 


(0.45) 


(4.57) 


(2.74) 


(1.69) 


(0.65) 


15000 


12.264 


8.999 


0.416 


5.317 


3.266 


1.833 


0.699 




(12.2) 


(8.98) 


(0.41) 


(5.31) 


(3.26) 


(1.84) 


(0.7) 


20000 


13.689 


9.998 


0.385 


5.922 


3.691 


1.944 


0.732 




(13.7) 


(9.98) 


(0.38) 


(5.91) 


(3.68) 


(1.94) 


(0.73) 
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TABLE II. Results for the energy of vortex states with n = 1, 2 and 3 of Rb atoms confined in an anisotropic harmonic 
trap with A = y/8 and w°/27r = 220/Ao Hz. Energy is in units of ftaj° . Numbers in the brackets correspond to the results 
obtained in the Thomas-Fermi limit. 



N 


{E!/N) K=1 


(Ei/N) K= 2 


(Ei/N) K= ; 


10 3 


4.455 


5.274 


6.164 




(3.385) 


(4.057) 


(4.402) 


5 x 10 3 


6.579 


7.237 


7.993 




(5.845) 


(6.513) 


(7.219) 


10 4 


8.176 


8.773 


9.469 




(7.537) 


(8.146) 


(8.853) 


5 x 10 4 


14.455 


14.94 


15.511 




(13.951) 


(14.397) 


(14.977) 


10 5 


18.78 


19.227 


19.771 




(18.306) 


(18.684) 


(19.192) 


5 x 10 5 


35.097 


35.474 


36.072 




(34.627) 


(34.876) 


(35.227) 


10 6 


46.134 


46.485 


46.888 




(45.637) 


(45.842) 


(46.134) 


5 x 10 6 


87.416 


87.72 


88.061 




(86.764) 


(86.892) 


(87.078) 


10 7 


115.229 


115.518 


115.837 




(114.459) 


(114.562) 


(114.714) 
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TABLE III. Results for the ground state of 7 Li atoms confined in an anisotropic harmonic trap with Ao = 0.72 and 
w^/27r = 163 Hz. Chemical potential and energy are in units of hu)± and length is in units a±. Numbers in the brackets 
correspond to the results of Ref. M . 



N 


Mi 


(Ei/N) 


{Ei/N) kin 


(E 1 /N) HO 


{E x /N) pot 


V< x l > 


V< 4 > 


1 


1.36 


1.36 


0.68 


0.68 





0.707 


0.833 






(1.36) 












100 


1.327 


1.344 


0.693 


0.670 


-0.017 


0.701 


0.824 


200 


1.291 


1.326 


0.707 


0.654 


-0.035 


0.695 


0.813 


300 


1.254 


1.309 


0.722 


0.641 


-0.054 


0.688 


0.803 


400 


1.214 


1.29 


0.74 


0.626 


-0.076 


0.681 


0.791 


500 


1.173 


1.271 


0.758 


0.611 


-0.098 


0.672 


0.786 


600 


1.125 


1.25 


0.782 


0.594 


-0.125 


0.665 


0.765 


700 


1.074 


1.229 


0.808 


0.576 


-0.155 


0.656 


0.75 


800 


1.017 


1.206 


0.839 


0.556 


-0.189 


0.645 


0.734 


900 


0.952 


1.182 


0.878 


0.533 


-0.23 


0.633 


0.715 


1000 


0.874 


1.155 


0.928 


0.507 


-0.28 


0.619 


0.693 






(1.15) 








(0.62) 


(0.69) 


1100 


0.776 


1.125 


0.999 


0.475 


-0.349 


0.60 


0.665 


1200 


0.625 


1.09 


1.121 


0.43 


-0.461 


0.573 


0.625 


1270 


0.346 


1.06 


1.42 


0.352 


-0.713 


0.521 


0.554 
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Figure Captions 



FIG. 1. Ground state energy per atom for 87 Rb as a function of N. Energy is in the units of TkUj_. The solid line is the 
result of our variational calculation. The dashed line is the result obtained by using the Gaussian trial wave 
function of Ref. |pd| , while the dotted line is the result obtained by using the Thomas-Fermi approximation. 
The filled circles are the results of Ref. || . 

FIG. 2. Ground state wave function for 87 Rb along the x axis for different values of N. The solid line is the result of 
our variational calculation. The dotted line is the result obtained by using the Gaussian trial wave function of 
Ref. jy]], while the dashed line is the result obtained by using the Thomas-Fermi approximation, (a) N — 100. 
(b)N = 1000. (c)N = 50000. 

FIG. 3. Aspect ratio in 87 Rb as a function of N. The lower and upper horizontal lines correspond to y/Xo and Ao 
respectively. 

FIG. 4. Vortex-state wave function of 5000 87 Rb atoms along the x axis for k = 1. 

FIG. 5. Ground state wave function for 7 Li along the x axis for N — 500 (lower curves) and N = 1270 (upper curves). 
The solid lines are the results of our variational calculation. The dotted line are the results obtained by using 
the Gaussian trial wave function of Ref. jll| . 

FIG. 6. Aspect ratio in 7 Li as a function of N. 
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